Interlayer hopping properties of electrons in layered metals 
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A formalism is proposed to study the electron tunneling between extended states, based on the 
spin-boson Hamiltonian previously used in two-level systems. It is applied to analyze the out-of- 
plane tunneling in layered metals considering different models. By studying the effects of in-plane 
interactions on the interlayer tunneling of electrons near the Fermi level, we establish the relation 
between departure from Fermi liquid behavior driven by electron correlations inside the layer and 
the out of plane coherence. Response functions, directly comparable with experimental data are 
obtained. 
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I. INTRODUCTION 

Layered materials have been object of intensive study 
since they present important physics. A key feature 
of many layered materials is the anisotropy exhibited 
by its transport properties: while being metallic within 
the layers, the transport in the c-axis, perpendicular to 
the layers, may be coherent or incoherent and undergo 
a crossover with temperature from one regime to the 
other, thus changing the effective dimensionality of the 
systemic. Correlation effects increase as dimensionality 
decreases, therefore dimensionality is crucial for the elec- 
tronic properties and to choose the appropriate model to 
study the system. Unusual properties are derived from 
the anisotropy and periodicity along the axis perpendicu- 
lar to the planes i.e. the structure of collective excitations 
absent in two dimensional (2D) and three dimensional 
(3D) electron gases^. 

Among the most studied layered materials are the 
high-temperature cuprate superconductors. These com- 
pounds present a strong anisotropy and are treated as 
two-dimensional systems in many approaches. In the 
normal state the transport properties within the doped 
Cu02 planes are very different from those along the o 
axis: electron motion in the c-direction is incoherent in 
contrast with the metallic behavior of the in-plane elec- 
trons as probed by the different p c and p a b resistivities 
and their different dependence with temperature^. Op- 
tical conductivity measurements confirm the anomalous 
c-axis properties^. These features have led to the asump- 
tion that the relevant physics of these substances lies on 
the CuC"2 planes, common to all the families. However, 
the importance of the c-axis structure is being intensively 
investigated since it offers a key to explain the differ- 
ences of the critical temperature T c and optimal doping 
in the different compounds. The relevance of the na- 
ture of the conductance in the direction perpendicular 
to the Cu02 planes to the nature of the superconduct- 
ing phase has been remarked on both theoretical^*^ and 
experimental 11 grounds. The anomalous behavior of the 
out of plane properties in the cuprates, and in analogy 



with the one dimensional Luttinger liquid, has led to the 
proposal of the failure of the conventional Fermi-liquid 
(FL) theory in these compounds^. 

Recent c-axis transport measurement focused on the 
resistivity anisotropy p c /p a b ratio indicate that p c /pa.b 
quickly increases with decreasing temperature. It was 
remarkably found that, at moderate temperature (be- 
tween 100 to 300K), p c / p a h is almost completely inde- 
pendent of doping x in the non superconducting regime 
(0.01 < x < 0.05), which suggest that the same charge 
confinement mechanism that renormalizes the c-axis hop- 
ping rate survives down to x = 0.01. On the contrary the 
Pc/pab ratio quickly changes when x enters into the su- 
perconducting regime, indicating that interlayer hopping 
quickly becomes easier as x is increased above 0.05, i.e. 
the charges are increasingly less confined in this doping 
ranged. Measurements of the out-of-plane magnetoresis- 
tance(MR) of La2_ x SR x CUO4 over a wide doping range 
and comparison with other systems suggest that a com- 
petition between the c-axis hopping rate and the in-plane 
scattering rate determines the behavior of the c-axis MR 
in layered perovskites with incoherent c-axis resistivity^. 

Many alternative models have been suggested to ex- 
plain the puzzling properties, from Fermi liquid modi- 
fied by strong electron-electron correlations^ to the non- 
Fermi Luttinger liquidt An alternative explanation of 
the emergence of incoherent behavior in the out of plane 
direction has been proposed in terms of the coupling of 
the interlayer electronic motion to charge excitations of 
the systemSiifl. This approach implicitly assumes that 
electron-electron interactions modify the in-plane elec- 
tron propagators in a non trivial way, at least at distances 
shorter than the elastic mean free path. The strong an- 
gular dependence of both the scattering rate and the in- 
terplane tunneling element can also lead to the observed 
anisotropics^. 

Graphite, another layered material, presents an intra- 
plane hopping much larger than the interplanc hybridiza- 
tion. Many of the transport properties established in 
the past for this well known material are being ques- 
tioned at present. Recent conductivity measurements re- 
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veal a suppresion of the c-axis conductivity much larger 
than what would be predicted by the band calculations 
of the interlayer hoppingi 6 -. Band structure calculations 
are also challenged by very recent claims of observation 
of quantum Hall plateaus in pure graphiteit The un- 
conventional transport properties of graphite such as the 
linear increase with energy of the inverse lifetime^ (see 
alsoi 9 -), suggest deviations from the conventional Fermi 
liquid behavior, which could be due to strong Coulomb 
interactions unscreened because of the lack of states at 
the Fermi leve l 20 ' 21 . 

By assuming that electron correlations modify the in- 
plane electron propagators^, we show that even in the 
clean limit, many-body effects can suppress the coherent 
contribution to the out of plane electron hopping. The 
clean limit is defined as that in which the elastic mean 
free path diverges. It is shown that for certain models 
of correlated electrons^Si* 2 ^ the interplane hopping be- 
tween extended states can be a relevant or an irrelevant 
variable, in the Renormalization Group (RG) sense, de- 
pending on the strength of the coupling constant. The 
scheme used here is based on the RG analysis as applied 
to models of interacting electrons 2 ^ 2 ^. We will demon- 
strate the usefulness of the method, by rederiving known 
results for a non trivial system, the array of Luttinger 
HquidaSSiSSiSi, and then we will apply the method to two 
well defined two-dimensional systems where ordinary per- 
turbation theory fails, due to the existence of logaritmic 
divergences. 

Recently there has been renewed interest in two- 
dimensional systems that support low-energy excitations 
which can be described by Dirac fermions. Examples 
are the flux-phase of planar magnet o 28 ! 29 , nodal quasi- 
particles in d-wave superconductor a 30 ' 31 , the insulating 
spin density wave phase of high T c superconductors 32 , 
and quasiparticles of planar zero-gap semiconductors as 
graphit e 20 ' 21 . Long range Coulomb interactions and dis- 
order are modelled in these systems by means of gauge 
fields coupled to the Dirac fermions. Most of these 
systems show anomalous transport properties ranging 
from mild departures of Fermi liquid behavior as in pure 
graphite, to the total destruction of the quasiparticle 
pole. Disorder can be modelled as random gauge fields 
coupled to the Dirac quasiparticles. It affects the quasi- 
particle Green's function in a computable way what in 
turn modifies the interlayer tunneling. 

This paper is organized as follows. In section II we 
present the method of calculation and show the results 
it gives in one dimension. In section III we apply it to 
two-dimensional models which show deviations from FL 
behavior: in particular to systems with a vanishing den- 
sity of states at the Fermi level as is the case of graphene 
sheets and to a system of planar electrons near a Van 
Hove singularity. The effect of disorder on these graphitic 
and related systems is discussed. The main physical con- 
sequences of our calculation are presented in section IV. 



II. THE METHOD OF CALCULATION. 

In the presence of electron-electron interactions, tun- 
neling processes are modified by inelastic scattering 
events. The influence of inelastic scattering on electron 
tunneling has been studied, using equivalent methods, 
in mesoscopic devices which show Coulomb blockade 33 , 
Luttinger liquida 25 i 26 i 27 i 34 , and dirty metals^. The sim- 
plest formulation of the method replaces the excitations 
of the system (such as electron-hole pairs) by a bath of 
harmonic oscillators with the same excitation spectrum. 
This approach can be justified rigorously in one dimen- 
sion, and is always an accurate description of the response 
of the system when the coupling of the quasiparticles to 
each individual excitation is weak, although the net ef- 
fect of the environment on the system under study can 
be large 3 ^. The expression for the coupling between the 
electrons and the oscillators is obtained by assuming that 
the oscillators describe the charge oscillations of the sys- 
tem. Then, the coupling can be related, using pertur- 
bation theory, to the charge-charge correlations in the 
electron gas. This is consistent with the assumption that 
the modes in the environment are weakly perturbed by 
their interaction to the low energy electrons whose tun- 
neling properties are been considered: 

Hint = etc; £^ (£)/>£ (1) 
k 

where cj(cj) creates(destroys) an electron at site i, and 
/5g describes the charge fluctuations of the environment, 
which are to be described as a set of harmonic modes as 
stated above. The interaction in Eq.Q is the simplest, 
and most common, coupling between the tunneling elec- 
trons and the excitations of the system, which is spin in- 
dependent. Other spin-dependent or more complicated 
couplings can also be taken into account, provided that 
the appropriate response function is used. 

Since the excitations of the system (electron- hole pairs, 
plasmons) are modelled as bosonic modes, one can write 
an effective electron-boson hamiltonian of the type: 

H-e-b = lielec + H-env + T~tint 

= ty^ta + u k b\b k + Y 9k,i4 c i( b t + W 

where TL e iec describes the individual quasiparticles, 7i en v 
stands for the set of harmonic oscillators which de- 
scribe the environment, and Hint defines the (linear) 
coupling between the two. The b\.{bk) are boson cre- 
ation(destruction) operators, the describe the elec- 
tronic hopping processes. The information about the 
interaction between the electron in state i and the en- 
vironment is encoded in the spectral function 36 

Ji(w) = Y lffc,i| 2( 5( w - ( 3 ) 
k 

The function Ji(ui) describes the retarded interaction in- 
duced by the environment. In standard perturbative 
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FIG. 1: Diagram corresponding to the lowest order contribu- 
tion to the self energy of an electron at site i in the hamilto- 
nian given in Eq.(|5J. This diagram can be used to describe 
the spectral function in Eq.lJSJ. 



treatments of the hamiltonian in Eq. J2Jl, this function is 
simply the self energy of an electron at site i, as schemat- 
ically shown in Fig. [p. 

Using second order perturbation theory and Eq. © , we 
can writo 33 ! 36 : 



J 4 H = 5>/(k)Im X (k, W ) 



(4) 



where x(k, oS) is the Fourier transform of the density- 
density response of the system, (pg(£)p_g(0)). The 
electron-boson interaction leads to a Franck-Condon fac- 
tor which reduces the effective tunneling rate. The 
Franck-Condon factor depends exponentially on the cou- 
pling between the particles and the oscillators, and when 
it diverges additional self consistency requirements be- 
tween the hopping amplitudes and the oscillator frequen- 
cies included in the calculation have to be imposed. The 
electron propagators acquire an anomalous time, or en- 
ergy, dependence, that can be calculated to all orders if 
the state i is localized, or, which is equivalent, neglecting 
the hopping terms in Eq.©. The present scheme can 
be considered a generalization to layered systems of the 
approach presented in reference^ for a Luttinger liquid, 
where the scaling of tunneling is expressed in terms of 
an effective density of states. The tunneling will become 
irrelevant, as inSi, in one of the models considered here, 
that of electrons in the vicinity of a saddle point of the 
dispersion relation, which is consistent with the picture 
of the c-axis interlayer tunneling strongly suppressed by 
voltage fluctuations, proposed by Turlakov and Legget^. 

Using the boson representation of the environment, we 
can write the electron propagator as: 



(4(t) Ci (f)) ~ (4(t) Ci (t')) x 



{-. 


j did 


I _ e M*-*') 


Ji(u)\ 








u> 2 J 



(5) 



exp • 



where (c|(t)cj(t'))o ~ e le ^* _ * is the Green's function in 
the absence of the interaction. This result can be viewed 
as the exponentiation of the second order self energy di- 
agram shown in Fig.© The approximation of the exci- 
tations in the environment by bosonic modes allows us 
to sum an infinite set of diagrams like the one in Fig. [Q . 



In contrast to previous work, we analyze tunneling be- 
tween coherent extended states. In order to do so, we 
need to generalize Eq.© to this case. We then need 
to know the Green's function of coherent states in the 
individual layers, G(k, w), including the correction due 
to the interaction to the environment. We firstly as- 
sume that Eq.© also holds in a system with extended 
states. For a standard metallic system, we must in- 
sert (cl(t)ci(t')) ~ l/(t - t') in Eq.©. It can be 
shown that this approximation is exact at short times, 
W <C (t — t')~ x <C A, where W is an energy scale related 
to the dynamics of the electrons, and A is the upper cutoff 
in the spectrum of the environment. This expression can 
be generalized, taking into account the spatial structure 
of the coupling to: 



<*1(W)> 



1 



t - 1 

x exp 



d 2 rdWkV C ( f - F '> 



k Jn Jn 



du> 



1-e' 



j(t-t') 



(6) 



where is the region of overlap of the wavefunctions 
Vf^r) and \&j(r). This expression, which can be seen as 
the exponential of the leading frequency dependent self- 
energy correction to the electron propagator, has been 
extensively used in studies of tunneling in zero dimen- 
sional systems (single electron transistors) which show 
Coulomb blockade^ 3 ., one dimensional conductors^, and 
disordered systems in arbitrary dimensions 3 ^. 

The effective interaction can, in turn, be written in 
terms of the response function as: 



U e// (k,^HU 2 (k)Im X (k, 



(7) 



The time dependence in Eq. © is determined by the 
low energy limit of the response function. In a Fermi 
liquid, we have 

Imx(k,u;) s» a(k)|a>| oj <C Ep , (8) 
where Ep is the Fermi energy. Eq. © then gives: 

lim (¥(t)$(t')) ~ l —— (9) 



where 



a= [ [ [ d 2 fd 2 f'rf 2 ke lC(F - ? 'V 2 (k)a(k) (10) 

The parameter a gives the correction to the scaling prop- 
erties of the Green's functions. Integration in k is re- 
stricted to |k| <c L^ 1 where L is the scale of the region 
where the tunneling process takes place. The value of L 
is limited by the length over which the phase of the elec- 
tronic wavefunctions within the layers is well defined. 



4 



The former approach is directly related to the theo- 
retical model proposed irA where the limit of coherent 
tunneling is set by the value a = 1. In sec. IIIB this pro- 
cedure is applyed to the study of interlaycr tunneling in 
layered systems whose Fermi surface is close to a saddle 
point which has been proposed as a possible explanation 
of some properties of the cuprates. 

We can also use Eq. © to analyze the interlayer tunnel- 
ing by applying Renormalization Group methods. The 
simplest case where this procedure has been used is 
for the problem of an electron tunneling between two 
states, i and j, which has been intensively studied^LS. 
We integrate out the high energy bosons, with energies 
A — dA < cok < A and rescaled hopping terms are de- 
fined. As mentioned earlier, Eq.© is valid for this range 
of energies. The renormalization of the hoppings is such 
that the properties of the effective Hamiltonian at en- 
ergies lo <C A remain invariant. If the hoppings ty- are 
small, any physical quantity which depends on them can 
be expanded, using time dependent perturbation theory, 
in powers of: 

The integration of the high energy modes implies that the 
terms in Ea. (|ll|l are restricted to t < A -1 , or, alterna- 
tively, the time unit have to be rescaled^, r' = re dA / A , 
where r ~ A -1 . Using Eq.©, the condition of keep- 
ing the perturbation expansion in powers of the terms in 
Eq. (|ll|l invariant implies that: 



(2+2q) 



(12) 



which can also be used to define the scaling dimension of 
the hopping terms. The beta function of the hopping is 
then 



0(WA) 
dl 



(13) 



log(Ao/A), and Aq is the initial value of the 



successfully used to 
in different situations 



where I 
cutoff. 

This approach has been 
describe inelastic tunneling 

The analysis which leads to Eq. (|13|) can be generalized 
to study hopping between extended states, provided that 
we can estimate the long time behavior of the Green's 
function, 



G(M-o = (4(t) C£ (0). 



(14) 



Equation(2|) is related to Ea.l|T5|l. by: 



lim G(r 

-?'|->0 



= / cTkG(k,w) 



(16) 



where D is the spatial dimension. In the cases discussed 
below, the interaction is instantaneous in time, and the 
non interacting Green's function can be written as: 



1 f h z 
G (k,w) <x-r[-s- 
uj \ to 



(17) 



where z = 1,2 depending on the dispersion relation of 
the model. In the following, we assume that the inter- 
acting Green's function has the same scaling properties, 
with the factor cl> _1 replaced by uj~ s in Eci. (|17fl . where 
S depends on the interactions. This can be shown to be 
correct in perturbation theory to all orders, in the models 
studied below, because the self energy diagrams have a 
logarithmic dependence on lu (it is a well known fact for 
the Luttinger liquid). Then, using eqs. (JTHJl, (|TB|l and 
(|17fl . we obtain: 



G(k,cu)oc\uj\ a ~ D / z ^[^- 



(18) 



and F(u) is finite. Thus, from the knowledge of the real 
space Green's function, using Eq.©, we obtain a, which, 
in turn, determines the exponent a — D / z which charac- 
terizes G(k,w). Generically, we can write: 



G, 



loc,ext 



(19) 



where the subindices loc, ext stand for localized and 
extended wavef unctions. In terms of these exponents, 
we can generalize Eq. <|13[) to tunneling between general 
states to: 



d(t l ° c > ext /A) 



01 



Uloc,ext ^ 



(20) 



Before proceeding to calculations of 5i oc and 5 ext for 
various models, it is interesting to note that, in general, 
the response function of an electron gas in dimension 
D > 1 behaves as 



lim 

*0,|k|->0 



so that, from Eq. lfTU|) . 



lim a ~ L^ 1 



We assume that, in a translationally invariant system, 
there is no dependence on the position of the local orbital, 
i. This result implies that the frequency dependence of 
the Green's function, in a continuum description, can be 
written as: 



lim G(r — r'jw) 

r — r'l— *Q 



(15) 



Thus, for D > 1, the contribution of the inelastic pro- 
cesses to the renormalization of the tunneling vanishes 
for delocalized states, L — * oo. This result is consistent 
with the existence of Fermi liquid behavior above one 
dimensioniS. Note that anisotropic Fermi surfaces, with 
inflection points, can lead to marginal behavior above 
one dimensioniiia^. 
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It is easy to show that, in an isotropic Fermi liquid in 
D dimensions, lim^^oo a oc L 1 ^ , where L is the lin- 
ear dimension of the (localized) electronic wavefunctions 
^(r). This result is due to the dependence on k of the 
response functions which goes as 

Imx(k, uj) 



fc^lkl 

Thus, for D > 1, we recover coherent tunneling in the 
limit of delocalized wavefunctions. 

In one dimension, one can use the non interacting ex- 
pression for Imxo (k, uj) , to obtain: 



d 2 fdWkV c ( F - F '> / 



j(t-t') 



U\ f0 t-t'«L/v F 
E F J \ log[v F (t - t')/L] t - f » L/v F W 



where we have assumed a smooth short range interaction, 
parametrized by U. Hence, the Green's functions have a 
non trivial power dependence on time, even in the L — > oo 
limit, in agreement with well known results for Luttinger 
liquids 34 . In order to obtain the energy dependence of 
the effective tunneling between k states near the Fermi 
surface, one needs to perform an additional integration 
over dr. In general, near a scale invariant fixed point, 
w oc |k| z , and for a ID conductor one knows that z = 1. 
Hence, 



Im G(w, k F ) 



OC UJ 



-z+a 



The flow of the hopping terms under a Renormalization 
Group scaling of the cutoff is^i: 



3(t/A) 
dl 



—a localized hopping 
1 — a extended hopping 



(22) 



where t denotes a hopping term, between localized or 
extended states. In the latter case, the hopping becomes 
an irrelevant variable 25 for a > 1. 



III. LAYERED MODELS 



A. Vanishing density of states at the Fermi level. 



In order to apply the procedure outlined in the previ- 
ous section, one needs the Fourier transform of the inter- 
action, 

y e// (k) = e 2 /(e |k|), 

where e is the electronic charge, and £q is the dielectric 
constant, and the susceptibility of the electron gas. For a 
single graphene plane, this quantity has been computed 
irt 2 ^ and is: 



Xo(k,w) 



(23) 



These expressions need to be inserted in equations (0 
and iJOJl- 

For simplicity, we consider the expression in Ea. (|23|l . 
as it allows us to obtain analytical results. The imagi- 
nary part, Imxo(k, uj), is different from zero if uj > u_F|k|. 
We cut off the spatial integrals at a scale, L, of the order 
of the electronic wavefunctions involved in the tunneling. 
Performing the same computation as in the case of D=l 
we obtain an expression similar to that in Eq. I|21|l except 
that the prefactor (U/E F ) 2 is replaced by the squared 
of effective coupling constant of the model, e 4 / (eoV F ) 2 . 
Thus, also in the graphene model, the propagators ac- 
quire an anomalous dimension what was advocated in3& 
as pointing out to a departure of the model from Fermi 
liquid behavior. As in ID, the value of the exponent z 
which relates length and time scales is z ~ 1. The scaling 
of the hoppings now are: 



A 



g(t/A) 
<9A 



1 



localized hopping 
extended hopping 



(24) 



The extra constant in the first equation with respect to 
Eq. (1221) reflects the vanishing of the density of states 
at the Fermi level for two dimensional electrons with a 
Dirac dispersion relation. 

In graphite, the dimensionless coupling constant, 
e 2 /v F , is of order unity. Under renormalization, it flows 
towards zero^fi. Thus, despite the departure from Fermi 
liquid behavior, interplane tunneling is a relevant vari- 
able and a coherent out of plane transport in clean 
graphene samples should be observed. The former pic- 
ture can change in the presence of disorder that is known 
to change the anomalous dimension of the fields^ 3 -. 



The simplest two dimensional model for interacting 
electrons where it can be rigourously shown that the 
couplings acquire logarithmic corrections in perturbation 
theory is a system of Dirac fermions (e^ = UF|k|), with 
Coulomb, 1/ 1 r — f'|, interaction. This model can be used 
to describe isolated graphene planes 20-21 , and can help 
to understand the anomalous behavior of graphite ob- 
served in recent experiments 1 ^^. Nodal quasiparticles 
in d-wave superconductors are also describable in these 
terms. 



B. Saddle point in the density of states. 

The Fermi surface of most hole-doped cuprates is close 
to a saddle point of the dispersion relationM. The possi- 
ble relevance of this fact to the superconducting transi- 
tion as well as to the anomalous behavior of the normal 
state was put forward in the early times of the cuprates 
and gave rise to the so called Van Hove scenario^ 5 ". 

We shall apply the mechanism described in section II 
to study the interlayer hopping of two electronic systems 
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described by the t— t' Hubbard model and whose Fermi 
surface lies close to a Van Hove singularity. 

The t-t '-Hubbard model has the dispersion relation 

e(k) = — 2t [cos(fc;j,a) + cos(k y a)] 

— 2t' cos(k x a) cos(kya) , (25) 

This dispersion relation has two inequivalent saddle 
points at A (jr, 0) and B (0, 7r). The Van Hove model in 
its simplest formulation is obtained by assuming that for 
fillings such that the Fermi line lies close to the singular- 
ities, the majority of states participating in the interac- 
tions will come from regions in the vicinity of the saddle 
points. In the present context the model is further rein- 
forced by the fact that the momentum dependence of the 
out of plane hopping in the cuprates tends to suppress 
tunneling from points away from the singularities. Tay- 
lor expanding Eq. 125JI around the two points gives the 
effective relation 



£A,s(k) « T(t T 2t')A£ ± (t ± 2t')k 2 y 



(26) 



The dynamical exponent (|17|l in this case is z=2. The 
dispersion relation (|26|l allows to formulate a renormal- 
izable effective model based on the hamiltonian: 

«i.fc<A 

^ = e i,k c k,i,s C i,k,s 
i=A,B;k,s 

+ ^ ] u i,i';s,s'Cj_k.s Ci ' ,s' c jit ,s" c i'" ,k"' ,s"' 7 (27) 

where Ao is a high energy cutoff which sets the limit of 
validity of the effective description. The particle-hole 
susceptibility has been computed m~: 



Imx(k,w) 



1 



47re(k) 



■e(k) 



E (k) 



(28) 



where e(k) is the dispersion relation 

The long time dependence of the Green's function is 
determined by the low energy behavior of x : 



lim Im x(k, us) 



UJ 



e(k) 



Inserting this expression in eqs.jfjj) and (JJJ, we can see 

that, irrespective of the details of the interaction, in the 
presence of a Van Hove singularity the exponent a in the 
time dependence of the Green's function goes as: 

lim a oc log(L) , 

£7 — >oo 

where L, as before, is the length scale which characterizes 
the wavefunction of the tunneling electron, and fl oc L 2 is 
the size of the integration region in Eq.©. The details of 
the anomalous dimension of the propagator in this case 
depend on the nature of the interactions which determine 
V eff (k,uj) = V 2 (k)Im X (k, uj) . 

To make contact with the work of Turlakov and 
LeggettS about the c-axis interlayer tunneling in the 



cuprates, we have computed the parameter a for differ- 
ent possible interactions. Ir.^ several models, describing 
the voltage noise, are examined to estimate the parame- 
ter a in order to describe the c-axis transport behavior. 
Here we consider three examples of potential: a long- 
range Coulomb interaction, a short-range Coulomb in- 
teraction and the intermediate situation between them, 
the Thomas-Fermi screened potential. The parameter a 
is calculated in the three cases. 

a) Unscreened Coulomb potential, VCk) — 2jL ~*-- Due 

to the highly singular interaction, the uj dependence of 
the effective potential is not linear and a is not well de- 
fined in this case. We can still analyze the tunneling by 
computing the effective potential Eq. 0. It is computed 
to be: 



fcgl 



1 + log 



uj 



where ujq = \Ep — Evh\ is a low-energy cutoff, Ep is 
the Fermi energy and the position of the saddle point is 
at Evh- keeps the chemical potential away from the 
singularity and is required to avoid infrared divergences 
in the integrals^. From Eq. I© it can be seen that the 
tunneling in this case is exponentially suppressed. 

b) Consider now a screened interaction of the Hubbard 
type, V(k) — Ua 2 , where a is the lattice unit. In this 
case the effective potential has a linear term in ui and 
the formalism described proceeds strightforward to give 
a value of the exponent a which is 



W 2 1 
a = ——am 



l V uj 



irt 2 (2tt) 2 
where A is a high-energy cutoff, and 



k x -\- ky 



Km = ^ log 



where k x ,k y ~ ^/uJo . The RG analysis ofiSdi shows that 
the coupling constant U of the model renormalizes to 
large values making the interlayer tunneling irrelevant. 

c) Finally we will analyze the case of a Thomas-Fermi 
screened potential which is an intermediate situation 
with respect to the previous cases. 



(2ne 2 ) 



2\2 



4(k 2 



k^p) 



In this case a can also be defined and it is computed to 
be 

?p 4 k 2 

— ■ ekTF \og 2 (k TF /k), 



e irt J 



where Utf is the Thomas-Fermi wavevector, and, k is 
a momentum cutoff k ~ y/u>o. In this case the longer 
range of the interaction makes the divergence softer as 
the dependence of the cutoff is as a squared log instead 
of a linear log. 
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1.75 
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FIG. 2: Effective potential calculated as a function of the 
energy for fixed k for a system of Van Hove layers coupled by 
Coulomb interaction . 

In all the three cases studied a diverges as Ep — > Evh- 
Thus, interlayer hopping is an irrelevant variable, and 
scales towards zero as the temperature or frequency is 
decreased. The additional logarithmic dependence found 
can be seen as a manifestation of the log 2 divergences 
which arise in the treatment of this model 22 . Note that, 
as in the graphene case, the coupling constants are also 
energy dependent, but in the Van Hove case they have an 
unstable flow and grow at low energies, suppressing even 
further the interlayer tunneling. This behavior of the 
interlayer hopping is in agreement with that obtained in 
Ref£ where it is found that c-axis interlayer tunneling is 
suppressed by voltage fluctuations. 

In the case of a layered system, we can use the Random 
Phase Approximation (RPA), and include the effects of 
interplane screening^: 

sinh(|k|d) 

^ '[cosh(|k|d) + ^ S inh(|k|d) X0 (k^)] 2 -1 

(29) 

where d is the interplane spacing. 

A numerical computation of 1291 provides the effective 
potential Eq. (JJJ. The imaginary part Xrpa gives the 
quantity usually known as the loss function that can be 
experimentally determined. This calculation is shown in 
Fig. 0. 

The calculated effective potential compares well with 
the experimental plots of the energy loss function of 
Bi2Sr2CaCu02, measured by transmission energy loss 
spectroscopy in the low energy region 4 ^, what reveals 
that the Van Hove model is also compatible with trans- 



port experiments. The results for the parameter a ex- 
tracted from the numerical computation are in qualita- 
tive agreement with the analytical expressions given in 
cases of b) screened Hubbard interaction and c) Thomas- 
Fermi screened potential, studied above. 

In a recent publication 48 it is argued that Umklapp 
processes are crucial to determine the spectrum of den- 
sity fluctuations, especially in two dimensional systems. 
It is worth to notice that the dispersion relation of the 
Van Hove model 12ri|) is a case where Umklapp processes 
enhance significantly the response function at low fre- 
quencies and wavevectors. 



IV. CONCLUSIONS. 

We have discussed the supression of interlayer tunnel- 
ing by inelastic processes in two dimensional systems in 
the clean limit. Our results suggest that, when pertur- 
bation theory for the in-plane interactions leads to loga- 
rithmic divergences, the out of plane tunneling acquires 
a non trivial energy dependence. This anomalous scal- 
ing of the interlayer hopping can make it irrelevant, at 
low energies, if the in-plane interactions are sufficiently 
strong. A well known problem where this non trivial scal- 
ing takes place is the tunneling between one dimensional 
Luttinger liquido 25 i 26 i 27 

In two dimensions, the scaling towards zero of the out 
of plane hopping is always the case if the Fermi level of 
the interacting electrons lies at a van Hove singularity 
(note that the Fermi level can, in certain circumstances, 
be pinned to the singularity 22 ). In this situation, in- 
sulating behavior in the out of plane direction is not 
incompatible with gapless or even superconducting in- 
plane properties. If the Fermi level is not tuned to the 
singularity, the scaling presented here is only valid for en- 
ergies larger than the distance of the Fermi level to the 
Van Hove singluarity. Within this range of energies or 
temperatures, insulating behavior can be expected, and 
Fermi liquid (coherent) behavior will set in at lower ener- 
gies. The appearance of instabilities can render this low 
energy regime unobservable. 

Clean graphene planes show the opposite behavior, as 
electron-electron interactions become irrelevant at low 
energies. 
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